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HYPERGEOMETRIC POLYNOMIALS ARE OPTIMAL 


D.V. BOGDANOV AND T.M. SADYKOV 


Abstract. With any integer convex polytope P C we associate a multivariate hy¬ 
pergeometric polynomial whose set of exponents is D P. This polynomial is defined 
uniquely up to a constant multiple and satisfies a holonomic system of partial differen¬ 
tial equations of Horn’s type. We prove that under certain nondegeneracy conditions the 
zero locus of any such polynomial is optimal in the sense of [7], i-e., that the topology 
of its amoeba m is as complex as it could possibly be. Using this, we derive optimal 
properties of several classical families of multivariate hypergeometric polynomials. 


1. Introduction 


Zeros of hypergeometric functions are known to exhibit highly complicated behavior. 
The univariate case has been extensively studied both classically (see e.g. mw and 
recently (see [3l[4l[2T] and the references therein). Already the distribution of zeros of poly¬ 
nomial instances of the simplest non-elementary hypergeometric function 2 F 1 (a, 6; c; x) is 
far from being clear. By letting the parameters a, 6, c assume values in various ranges, 
one can obtain a wide variety of shapes. Some of them are highly regular (see e.g. Fig. 
while other are nearly chaotic. 



Figure 1. The hypergeometric aster: zeros of the polynomials 



Polynomial instances of hypergeometric functions in one and several variables are very 
diverse. They comprise the classical Chebyshev polynomials of the first and the second 
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kind, the Gegenbauer, Hermite, Jacobi, Laguerre and Legendre polynomials as well as 
their numerous multivariate analogues j5]. 

Despite the diversity of families of hypergeometric polynomials, most of them share the 
following key properties that justify the usage of the term ” hypergeometric”: 

1. The polynomials are dense (possibly after a suitable monomial change of variables). 

2. The coefficients of a hypergeometric polynomial are related through some recursion 
with polynomial coefficients. 

3. For univariate polynomials, there is typically a single representative (up to a suitable 
normalization) of a given degree within a family of hypergeometric polynomials. 

4. All polynomials in the family satisfy a differential equation of a fixed order with 
polynomial coefficients (or a system of such equations) whose parameters encode the 
degree of a polynomial. 

5. In the case of one dimension, the absolute values of the roots of a classical hy¬ 
pergeometric polynomial are all different (possibly after a suitable monomial change of 
variables). 

6. Many of hypergeometric polynomials enjoy various extremal properties. 

In the present paper, we introduce a definition of a multivariate hypergeometric poly¬ 
nomial in n > 2 complex variables that is coherent with the properties 1-6 listed above. 
Namely, with any integer convex poly tope P C M” we associate a multivariate hyper¬ 
geometric polynomial whose set of exponents is n P. For this polynomial to be truly 
hypergeometric in the sense made precise below, we need to assume that any pair of points 
in Z"nP can be connected by a polygonal fine with unit sides and integer vertices. This 
assumption does not affect the generality of the results since any polytope that does not 
satisfy this condition gives rise to a finite number of hypergeometric polynomials that 
can be considered independently. The assumption of convexity of the polytope P is then 


automatically satisfied as clarified in Lemma 3.4 


The hypergeometric polynomial associated with the polytope P is defined uniquely up 
to a constant multiple and satisfies a holonomic system of partial differential equations of 
Horn’s type [161 [H]. We prove that under certain nondegeneracy conditions (see Theo- 
3.10) the zero locus of any such polynomial is optimal in the sense of [7|. Generally 


rem 


speaking, this means that the topology of the amoeba of such a polynomial is as 

complicated as it could possibly be (see Definition |2.8| ). This property is the multivariate 
counterpart of the property of having different absolute values of the roots for a poly¬ 
nomial in a single variable. We show various families of classically known multivariate 
polynomials to be optimal: a biorthogonal basis in the unit ball, certain polynomial in¬ 
stances of the Appel Fi function, bivariate Chebyshev polynomials of the second kind etc. 
We also introduce a simplicial complex that encodes intrinsic combinatorial properties of 
an algebraic variety while possessing key properties of its compactified amoeba. 

Pictures of amoebas in the paper have been created with Matlab 7.9. The authors are 
thankful to L. Lang for valuable comments on the paper. 
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2. Hypergeometric systems and AMOEBAS 


Throughout the paper, we denote by n the number of x G C** variables. For a = 
(«!,..., an), we use the notation |q;| = a\ = ai \... «„!. For x = (xi,..., x„) 

and a = {ai ,..., «„), we denote by x" the monomial x“^ ... x"”". 

Definition 2.1. A formal Laurent series 

( 2 . 1 ) 

sez" 

is called hypergeometric if for any j = 1,... ,n the quotient (p{s + ej)/(p{s) is a rational 
function in s = (si,..., s^). Throughout the paper we denote this rational function by 
Pj{s)/Qj{s + Cj). Here is the standard basis of the lattice Z”. By the support of 

this series we mean the subset of Z"" on which (y9(s) jtz Q. 

By a hypergeometrie funetion we will mean a (typically multi-valued) analytic function 
obtained by means of analytic continuation of a hypergeometric series with a nonempty 
domain of convergence along all possible paths in C**. 


Theorem 2.2. (Ore, Sato, see [8].) The coefficients of a hypergeometric series are given 
by the formula 

m 

(2.2) ip{s)=fU{s) J]r((A,,s} + Q), 

1=1 

where P — tff ... t^, U, Ci G C, Aj = (Aj^i,..., Aj^„) G Z**, i = 1,..., m, and U(s) is the 
product of a certain rational function and a periodic function (^{s) such that 4>{s + Cj) = 
(j){s) for every j = 1,... ,n. 


Given the above data (tj, q, Aj, {7(s)) that determines the coefficient of a hypergeomet¬ 
ric series, it is straightforward to compute the rational functions Pi{s)/Qi{s + ef) using 
the F-function identity. The converse requires solving a system of difference equations 
which is only solvable under some compatibility conditions on Pi,Qi. A careful analysis 
of this system of difference equations has been performed in [16] . 


We will call any function of the form (2.2) the Ore-Sato coefficient of a hypergeometric 
series. In this paper the Ore-Sato coefficient (2.2) plays the role of a primary object 
which generates everything else: the series, the hypergeometric system of differential 
equations, its polynomial solution (if any) and its amoeba. We will also assume that 
m > n since otherwise the corresponding hypergeometric series (2.1) is just a linear 
combination of hypergeometric series in fewer variables (times an arbitrary function in 
remaining variables that makes the system non-holonomic) and n can be reduced to meet 
the inequality. 


Definition 2.3. The Horn system of an Ore-Sato coefficient. A (formal) Laurent series 
SseZ" whose coefficient satisfies the relations ip{s-\-ej)/ip{s) = Pj{s)/Qj{s-\-ej) is a 

(formal) solution to the following system of partial differential equations of hypergeometric 
type 


(2.3) 


= QiV)f(^). i = 1..... 


n. 
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Here 9 = (6*i,..., On), Oj = The system (2.3) will be referred to as the Horn hyperge¬ 

ometric system defined by the Ore-Sato coefficient ip{s) (see [H]) and denoted by Horn((/7). 
In this paper we only treat holonomic Horn hypergeometric systems, i.e. rank(Horn((y?)) 
is always assumed to be finite. 


We will often be dealing with the important special case of an Ore-Sato coefficient (2.2) 
where U = 1 for any z = 1,..., n and U{s) = 1. The Horn system associated with such an 
Ore-Sato coefficient will be denoted by Horn(74,c), where A is the matrix with the rows 
Ai,..., Am G 'Zfi and c = (ci,..., Cm) G C"*. In this case the following operators Pj{0) 
and Qj{9) explicitly determine the system (2.3): 

n n 


— n H ((A,, s)+c, + 








Aj, s) -I- Cj -I- 




Definition 2.4. The support of a series solution to (2.3) is called irreducible if there 
exists no series solution to (2.3) supported in its proper nonempty subset. 

Definition 2.5. The amoeba Af of a Laurent polynomial f{x) (or of the algebraic hy¬ 
persurface {f{x) = 0}) is defined to be the image of the hypersurface /“^(O) under the 
map Log : {xi,...,Xn) ^ (log |a:i|,..., log |x^|). 


Despite losing n real dimensions, the amoeba of an algebraic hypersurface encodes its 
several intrinsic properties [zmann]. The main results of the paper describe topological 
properties of the amoebas of hypergeometric polynomials. The next lemma shows that 
certain transformations of a polynomial do not affect the topology of its amoeba. 


Lemma 2.6. The number of conneeted eomponents of the amoeba complement of a (Lau¬ 
rent) polynomial p{xi ,..., Xn) is the same as that of the polynomial x^'pftix "^, • • •, tnX'^'^Y 
for any £ G N, a = (ui, ..., a„) G 'IT' , t = (ti, ... ,tn) G (C*)” and any nondegenerate 
integer matrix v with the rows ui, ... ,n„. That is, there is a bijection between the con¬ 
nected components of the complements of the two amoebas; moreover, the orders [7j and 
the recession cones [13] of the corresponding components are transformed into each other 
by the linear map with the matrix v. 


Proof. A monomial factor can only vanish in the union of the coordinate hyperplanes 
that is mapped off to infinity by the logarithmic map. The amoeba does not reflect the 
multiplicities of the zeros of a polynomial and it can therefore be raised to any positive 
power. The map (xi,... ,Xn) (tiXi ,... ,tnXn) corresponds to the shift of the amoeba 
space with respect to the vector log |ti|,... ,log |t„|. Finally, a nondegenerate monomial 
change of variables in the complex torus (C*)” corresponds to the linear transformation 
of the amoeba space defined by the matrix of exponents of the monomials. Clearly a 
nondegenerate linear map preserves the topology of amoebas and provides a bijection for 
the recession cones of the complement components. The last statement of the lemma is 
an immediate consequence of the definition of the order of a component in the amoeba 
complement [7| . □ 

Recall that the Newton polytope Mp(x) of a Laurent polynomial p{x) is defined to be 
the convex hull in of the support of p{x). The following result shows that the Newton 







HYPERGEOMETRIC POLYNOMIALS ARE OPTIMAL 


5 


polytope Mp{x) reflects the structure of the amoeba Ap{x) [3 Theorem 2.8 and Proposi¬ 
tion 2.6]. 

Theorem 2.7. (See [7].) Letp{x) be a Laurent polynomial and let {M} denote the family 
of connected components of the amoeba complement ‘^Ap(x) ■ There exists an injective func¬ 
tion u : {M} ^ F n Ap(a;) such that the cone whieh is dual to Mp(x) the point u(M) 
coincides with the recession eone ra of M. In partieular, the number of eonnected eom- 
ponents of^Ap(x) cannot be smaller than the number of vertices of Np(x) o,nd cannot exceed 
the number of integer points in Mp{^x) • 

The two extreme values for the number of connected components of the complement of 
an amoeba are of particular interest. 

Definition 2.8. (Cf. [3, Dehnition 2.9].) An algebraic hypersurface H C (C*)”^, n > 2, 
is called optimal if the number of connected components of its amoeba complement ^An 
equals the number of integer points in the Newton polytope of the defining polynomial 
of H. We will say that a polynomial (as well as its amoeba) is optimal if its zero locus is 
an optimal algebraic hypersurface. 


Since the amoeba of a polynomial does not carry any information on the multiplicities 
of its roots, any one-dimensional amoeba (which is just a finite set of distinct points in 

k 

tti,... ,ak G M) can be treated as the amoeba of the polynomial H of whose 

i=i 

roots are positive and distinct. Thus Definition |2.8| is trivial in the univariate case. The 
correct extension of Definition |2.8| to one dimension is to say that all the roots of the 
polynomial in question have different absolute values. 

Example 2.9. In accordance with Definition |2.3| the Ore-Sato coefficient 

(fi{s, t) = r(—s — t i)r(2s — t — 2)r(—s -i- 2t — 2 ) 

yields the polynomials 

Pi{s, t) = (2s — t — 2)(2s — t — 1), Qi{s,t) = {s — 2t 2){s 1 — 1), 

-^ 2 ( 5 , t) = (s — 2t -|- l)(s — 2t -|- 2), Q 2 {s, t) = —(2s — t — 2)(^s 1 — 1). 

The corresponding Horn hypergeometric system is given by the linear differential operators 

T (2ex -Op- 2)(2ex -ey-i)- (6x - 2d. + 2)(d^ + -1), 

y { 9 x - 29 y + l)(d„ - 29 y + 2) + (2d„ - 9 y - 2)(d^ 9 - 9 y - 1). 

It is straightforward to check that ( 2.4[ ) is satisfied by the hypergeometric polynomial 
p{x, y) = X y 6xy x^y'^. 

The amoeba of p{x, y) together with its compactified version [H] are depicted in Fig.[^ 
We remark that the actual shapes of the amoeba and the compactified amoeba of p(x, y) 
(with respect to both directions of the amoeba tentacles and the complement of its com¬ 
pact version) are rather different from those given in |llj . 


(2.4) 


An algebraic hypersurface is optimal if the topology of its amoeba is as complicated 
as it could possibly be in the view of Theorem |2.7| (that is, the number of connected 
components in the amoeba complement is maximal). A two-dimensional optimal amoeba 
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Figure 2. The amoeba and the compactified amoeba m ofp{x,y) 


has the maximal possible number of bounded connected components in its complement 
and the maximal numbe r of parallel tentacles. 

we will not distinguish polynomials whose zero loci in (C )'^ 


In the view of Lemma 


2.6 


can be transformed into each other by a nondegenerate monomial change of variables. 
The reason for this is illustrated by the following example. 


Example 2.10. A typical example of a family of optimal polynomials arising in a different 
theory is given by the biorthogonal family in the unit ball {14(a:)}„gi^n , x G C” defined 
through their generating function (see [5l Section 2.3]): 

(1 — 2(a,x) + ||o|p)^“ = a“I4(x). 

aGNg- 


Neglecting an inessential monomial factor of 14 (x) (whose zero locus is contained in 
the uni on o f the coordinate hyperplanes and therefore does not affect the amoeba, see 
Lemma 2.6) one can represent Va{x) = 14(0 with 14 being a polynomial in 4 = x]- One 


re 


can check that the zero locus of 14(0 nn optimal algebraic hypersurface in C”. T 
Newton polygon and the amoeba of the bivariate polynomial 0(6,lo) are depicted in Fig. ^ 
Numerous other families of optimal multivariate polynomials can be found in [HI Cha 
ter 2]. 


3. Hypergeometric polynomials in several variables 


There exist several well-known and important families of hypergeometric polynomials 0 
El]. What they all have in common is the fact that they satisfy certain linear differential 
equations with polynomial coefficients that are special instances of (2.3). However, the 
family of all holonomic systems of partial differential equations of the form (2.3) is far too 
vast to serve as a definition of a hypergeometric polynomial. In fact, from the point of 
view of the general Definition |2.1| given above, any polynomial in any number of variables 
is a hypergeometric function. This is made precise in the following statement. 


Theorem 3.1. For any polynomial p{x) G C[xi,... ,xj there exists a noneonfluent [T3| 
holonomie hypergeometrie system of the form ( fA3| ) having p{x) as one of its solutions. 
Moreover, it ean he chosen in such a way that the hypergeometrie ideal in the Weyl algebra 
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Figure 3. The Newton polygon of V{ia,io){x^y) in [HI Section 2.3] and its amoeba 


defining this system admits a basis that consists of a commutative family of differential 
operators. 


Proof. In fact, any given polynomial is annihilated by a family of hypergeometric ideals 
that contains continuously many elements. To prove the theorem, we will present an 
explicit representative with desired properties. 

Let p{x) = ^ CaX°‘ be a polynomial with the support S. Here and throughout the proof 

aes 

we assume S to be finite. We denote by the cardinality of S and let |q;| =«! + .. . + «„• 
For s G define the Ore-Sato coefficient (/p(s) by 


(3.1) 




]]][ (si + . . . + — |tt|) 

ctGS' 

n 

n n i^j ~ ^j) 

j=l aeS 


By Definition |2.3[ the action of the j-th hypergeometric differential operator in the system 
defined by this Ore-Sato coefficient on p{x) is given by 



since C[0u...,0n] is a commutative subring in the Weyf algebra and (6*i + ... + — 

| q ;|) t “ = (6j — afix'^ = 0 for any j = 1 ,..., n. 

The hypergeometric system defined by the Ore-Sato coefficient ( |3.1| ) is nonconfluent 
by definition. By Theorem 2.8 in ini this system is hofonomic with the holonomic 
rank (fifSY. The fact that it is generated by a commutative famify of hypergeometric 
operators foflows from Lemma 2.5 in [m. □ 
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Since every monomial in p[x) is annihilated by each operator in the hypergeometric 
system defined by (3.1), the same argument works for any Puiseux polynomial with arbi¬ 
trary exponents in C”. Although p(x) is in the kernels of the differential operators that 
form a holonomic hypergeometric system, the monomials in p{x) are in no way related 
to each other. To obtain a meaningful definition of a hypergeometric polynomial based 


on (2.3), we will impose further assumptions on the Ore-Sato coefficient that defines the 


system. 

Throughout the rest of the paper we will only consider polynomials in n variables 
whose Newton polytopes have nonzero n-dimensional volume. For if the volume of such a 
Newton polytope is zero, a suitable monomial change of variables can be used to reduce 
the number of variables. 

From now on we will adopt the following definition. 

Definition 3.2. A set A C Z” is called 'IT'-convex if the condition + (1 “ : 

A G [0,1]} n IT C S holds for any G S. 

Definition 3.3. A set A C IT is said to be connected if any two points of this set 
can be connected by a polygonal line with unit sides and vertices in S. 

For insta nce, the support of the bivariate hypergeometric polynomial x-\-y-\-Qxy-\-x‘^y'^ 


in Example 


2.9 


IS a 


-convex but not a I -connected set. In fact, it consists of two 


-connected components: {(1, 0), (0,1), (1,1)} and {(2,2)}. 


Recall that the support A of a solution to the hypergeometric system (2.3) is called 
irreducible if there is no nonzero solution to ( 2.3[) sup ported in a proper nonempty subset 
of S. Any irreducible support of a solution to (2.3) is always a Z”-connected set. The 
following statement has been established in 


Lemma 3.4. (See [T7]. ) If the support S of a polynomial solution to the system is 
irreducible, then S is a -convex set. 


We next show that any convex integer polytope supports an irreducible solution to a 
suitable instance of the hypergeometric system (2.3). 

Lemma 3.5. For any convex integer polytope P G M” such that PflZ"" is IT-connected, 
there exists a hypergeometric system of the form (2.3) and its polynomial solution p(x) 
with irreducible support such that Mp(x) = P- 


Proof. Let {Bj, s) Cj = 0, j = 1,..., g be the equations of the hyperplanes containing 
the faces of P with Bj being the outer normal to P at the respective face. Since P is 
an integer polytope, we may without loss of generality assume the components of the 
vector Bj to be integer and relatively prime. 

Consider the Ore-Sato coefficient 

q 

p{s) = + 9 )- 

By Definition |2.3| the hypergeometric system defined by (p{s) only depends on the quotients 
Rj{s) := ip{s + ej)/ip{s) that are rational functions in s. Using the F-function identity 
F(z) = together with the fact that the meromorphic function is periodic 
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with the period 1, we conclude that the quotients Rj{s) coincide with those for the entire 
function 


(3.2) 


exp I ITT I +Cj 


The support of the polynomial with the coefficient (3.2) does not change if we replace 
the numerator of (3.2) by 1. In fact, this numerator is the exponential part of the 


Ore-Sato coefficient (2.2) and by Lemma 2.6 affects neither the support of the polynomial 


solutions to the corresponding hypergeometric system nor the topological properties of 
their amoebas. For these reasons we define the coefficient of the hypergeometric polyno¬ 
mial under construction to be 


(3.3) i/>p(s) := -. 

nr(l-(Si.«)-<;,) 

i=i 

The function 'ipp{s) is completely defined by the integer polytope P. By the construction 
the function 'ipp{s) vanishes at any lattice point that does not belong to P and is positive 
in P. Define the polynomial p{x) to be 

p{x) = tpp{s)x^. 

By the explicit construction, the polynomial p{x) is supported in P fl TT and satisfies 
the hypergeometric system Horn('^p(s)). The support is irreducible since P fl 'IT is Z”- 
connected and since 'ipp{s) does not vanish in P. □ 

The conclusion of the above lemma still holds even without the condition of 
connectedness of the set of integer points in the defining polynomial. However, the support 
of the polynomial produced by the construction in the proof of the lemma will no longer 
be irreducible. Such a polynomial cannot be considered as ’’truly hypergeometric” since 
it is a linear combination of two or more polynomials satisfying the same hypergeomet¬ 
ric system, see Example |2.9[ The properties of the amoeba of such a polynomial are in 
general heavily dependent on the coefficients of this linear combination. Yet, it is always 
possible to choose these coefficients in such a way that the hypergeometric polynomial 
with the (reducible) support P fl Z” is optimal. Thus we may and will without loss of 
generality assume throughout the rest of the paper that the set P fl IT' is Z’^-connected. 


Remark 3. 6. O ne can still not define a hypergeometric polynomial to be a polynomial 
solution to (2.3) with a Z’^-convex irreducible support S since any polynomial supported 
in S will satisfy this condition. This can be seen by introducing more factors into the 
Ore-Sato coefficient that will affect the coefficients of the polynomiai soiution but wifi not 
corrupt its support. Instead, we wili distinguish the oniy polynomiai that has support S 
and satisfies the hypergeometric system of the smaliest possibie hoionomic rank. 
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The following definition is central in the paper and brings together the intrinsic proper¬ 
ties of the classical families of hypergeometric polynomials: the denseness of the support, 
the irredncibility of the support and the property of being a solution to a suitable system 
of linear differential equations with polynomial coefficients. 


Definition 3.7. By a multivariate hypergeometric polynomial supported in a convex in¬ 
teger polytope P G n > 2 we will mean the polynomial 


(3.4) 




sePnZ” 


with 'ipp{s) defined by (3.3). 

By construction, a translation of the defining polytope P by an integer vector results in 
multiplication of the corresponding hypergeometric polynomial with a monomial. Thanks 
to Lemma 2.6 this does not affect the amoeba of (3.4). Throughout the rest of the paper 


we will identify polytopes that are translations of each other with respect to an integer 
vector. 


Remark 3.8. The hypergeometric polynomial introduced in Definition |3.7| satisfies the 
hypergeometric system of the smallest possible holonomic rank among all hypergeometric 
systems that admit an irreducible polynomial solution with the support P. This property 
can be used as the definition of a P-supported hypergeometric polynomial. 


In one dimension, Definition |3.7| yields a class of polynomials that is far too small 
to be interesting. Namely, for a segment [a, 6] C M with the inte ger e ndpoints a, b the 
polynomial x‘^{x + 1)^““ is hypergeometric in the sense of Definition 3.7, satisfies a hyper¬ 
geometric differential equation (that is a special instance of (2.3)) of the smallest possible 
holonomic rank 1 and vanishes on an optimal algebraic set. In what follows we will focus 
on the multivariate case. 


Example 3.9. Here we compute hypergeometric polynomials associated with certain 
families of integer convex polytopes and investigate their properties. 

1) Direct product of segments. If the polytope P in Definition 3.7| is the direct product of 
segments, we may without loss of generality assu me it to be P = [0, ai] x [0, a 2 ] x... x [0, a„] 


for Oj G N. In this case the Ore-Sato coefficient (3.3) is given by 


IKn 


^3 + - Sj + 1 )) 


-1 


i=i 


The corresponding hypergeometric polynomial is a constant multiple of U + 1)“^ . Its 

i=i 

amoeba is the union of the coordinate hyperplanes in M and is optimal in the sense of 
Definition 12.81 

2) A simplex. Let now the polytope P C M” be defined as the convex hull of the origin 
and the points (0,..., fc,..., 0) (fc in the j-th position) for j = 1,..., n. The corresponding 
Ore-Sato coefficient is given by 



A: -I- 1 - ^ Sj I ]q r{sj + 1) 


i=i 


i=i 


-1 
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its amoeba being just a hyperplane 


The hypergeometric polynomial defined by this co effic ient is a constant multiple of [xi + 
.. .+x„+l)^. It is optimal in the sense of Definition 21 
amoeba [7j. 

3) Cross-poly topes. Recall that the n-dimensional cross-polytope is the convex polytope 
whose vertices are all the permutations of (±1, 0,..., 0) G The only integer point of 
a cross-polytope that is not its vertex is the origin. The Ore-Sato coefficient defined by 
the n-dimensional cross-polytope is given by 


-1 


n 


r(6:i5i + . . . + £nSn + 2) 


, £j=zbl, j = 


It might happen that an algebraic hypersurface is defined by a polynomial which 
is irreducible in C[ti,. .. but can be factored in the ring of Puiseux polynomials 
for some d G N (see definition of strong irreducibility below). If both 
factors represent branches of the same hypersurface, such a polynomial cannot be op¬ 
timal. It turns out that a polynomial of this kind can be hypergeometric in the sense 
of Definition |3.7[ For instance, let n = 2 and the polytope P be the two-dimensional 
cross-polytope with the vertices (1, 0), (0,1), (—1,0), (0, —1). The Ore-Sato coefficient as¬ 
sociated with this polytope is 

(r(2 — s — t)r(2 — s -i- t)r(2 -i- s — t)r(2 -i- s -i- 1))~^ 


and the corresponding hypergeo metri c polynomial is given (up to a monomial multiple 
that is unimportant due to Lemma 2.6) by p(x, y) = x-by+Axy+x^y+xy"^. This polynomial 
is not optimal. This can be seen by the direct computation of its amoeba. An alternative 
way to prove this is to observe that the origin must belong to the component of order 
(1,1), if any, in the complement of the amoeba of the polynomial x-\-y + axy -\- x^y xy"^. 
However, p(—1, —1) = 0 while the point (—1, —1) is mapped to the origin by the map Log. 
In fact, p{x,y) is ”on the boundary” of the set of optimal polynomials supported in P. 

The above becomes clear in view of the following factorization of p{x, y) in Cl[^/x, ^Jy] : 
p{x, y) = (^/^(l -\-y) + \/^^/y{l + a:)) {\/x{l -\-y) - + a:)) . Both factors rep¬ 

resent branches of the same algebraic hypersurface {p(x, y) = 0} which can therefore not 
be optimal. 

n 

We remark that the generic polynomial p{xi, ..., Xn) = c-\- ^ 


i=i 


OjXj 


— ) supported in 


the n-dimensional cross-polytope and having positive coefficients Oj, bj, c is optimal if and 

n 

only if ^ \/(ijbj < |. Thus the hypergeometric polynomial defined by the n-dimensional 
i=i 

cross-polytope is optimal if and only if n > 2. 

4) The Hirzebruch surface. Recall that the Hirzebruch surface Fi is defined by the fan 
generated by (1, 0), (0,1), (—1,1) and (0, —1). The hypergeometric polynomial supported 
in (a translation of) the convex hull of these vectors is defined by the Ore-Sato coefficient 


(r (3 — t)r (4 — s — t)r(t — s -h 2)r(2s -i -1 — i)) 


-1 


This polynomial is a constant multiple of 3x -T I2xy -\- 2x^y 2y^ 3xy^ and is optimal. 
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(a) 


(b) 


(c) 


(d) 


Figure 4. The amoeba, compactified amoeba, the weighted compactified 
amoeba of the 6th Hadamard power and a vanishing connected component 


of the complement to the WCA(p 3 (x,y)) (see Definition 3.12) 


Recall that a (Laurent) polynomial p{xi,... ,Xn) E C[a:f ,... ,x^^] is called strongly 
irreducible if for any nondegenerate integer matrix v the polynomial p{x'^^, ■ ■ ■, x'"'^) is irre¬ 
ducible over C. Here vj is the jth row of the matrix v. This definition is only meaningful 
in n > 2 variables since there are no univariate strongly irreducible polynomials. The 
property of a polynomial being strongly irreducible in ..., x^^] is a generic condi¬ 

tion in the Zariski topology. Example 3.9| 3) provides a polynomial that is irreducible but 
is not strongly irreducible. 

By Lemma |2.6| the amoeba of a polynomial that is not strongly irreducible is the union 
of the amoebas of its Puiseux polynomial factors. We thus may without any loss of 
generality restrict our attention to strongly irreducible hypergeometric polynomials. 

The hypergeometric polynomials introduced in Definition |3.7| enjoy properties that are 
parallel to the properties 1-6 listed in the Introduction of the classical hypergeometric 
polynomials. Indeed, for any integer convex polytope there is only one (up to a constant 
multiple) hypergeometric polynomial of the form (3.4) supported in this polytope; this 


polynomial is dense; it satisfies the hypergeometric system Horn('i/ip(s)). The counterpart 
of the properties of the roots of (3.4) is given by the next theorem which is the main result 
of the paper. 

Theorem 3.10. Strongly irreducible hypergeometric polynomials are optimal. 


Proof. Let p{x) be a hypergeometric polynomial in the sense of Definition 3.7 with the 
variables x = {xi,... ,Xn) E C” and supported in a finite Z”^-connected set S C Z”. 
Denote by (p(s) the Ore-Sato coefficient of p{x). By definition, the function (p(s) is well- 
defined, finite and positive on S. 

By the comment after Remark |3.8| a univariate polynomial that is hypergeometric in 


the sense of Definition 3.7 is optimal. From now on we will consider the case n> 2. 

By0, the order of the connected component M of the complement of the amoeba Af 
is given by the vector with the coordinates 

1 f Ojf(x) dxi A ... A dxn 

' J = l,...,n. 


{2m) 


Log h€) 


fix) 


Xi .. .x„ 


Here ^ is any point in M. We will show that for any a E S there exists a con necte d 
component of order a E S C Z” in the amoeba complement b4p(a;). By Lemma 


2.6 


it 
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suffices to consider the case when a belongs to the interior of the convex hull (in of S. 


Indeed, if « is a vertex of Np^x) then the desired conclusion is established in Theorem 2.7 
If a is on the boundary of the convex hull of S but is not its vertex, then a suitable 
nondegenerate monomial change of variables can be used to reduce the dimension of the 
variable space and put a into the interior of the Newton polytope of a polynomial in fewer 
variables obtained from p{x) by setting some of its variables to zero. 


Furthermore, we may without loss of generality assume that a = 0. Indeed, by Lemma [2i6 
the amoebas of p{x) and x~^p{x) are the same and the component of order a in ^Ap[x) 
coincides with the component of order 0 in ^Ax-c ^pjx) - 

By the Bohr-Mollerup theorem the coefficient ( |3 3[ ) of a hypergeometric polynomial is a 
positive and strictly logarithmically concave function on the convex hull of the support of 
this polynomial. It follows from [T5l Corollary 2] that a sufficiently high positive (integer) 
power ^ of the Ore-Sato coefficient <p(s) defines an optimal hypergeometric polynomial. 
In fact, any sufficiently big positive real power of (p(s) will define an optimal polynomial 
but it will only satisfy a hypergeometric system of equations for integer powers. 

It remains to check that £ = 1 works. Let the Ore-Sato coefficient 'ipp{s) be defined 
by (3.3). By definition the function 




ses 
s / 0 


is positive and convex in M” . It tends to infinity as xj ^ 0 or ^ oo for any j = 1,... ,n. 
Thus there exists the unique minimum p = (pi,... ,Pn) G IR” such that p{x) > p(p) for 
any x G . 

Definition 3.11. Following the ideas of [20], we define the weighted moment map asso¬ 
ciated with the algebraic hypersurface {x G C” : /(t) := ^ = 0} through 

s&S 


Pf{x) 


J2 ^ 

ses 

ses 


It follows from the general theory of moment maps [0] that /i/(C”) C A//. 


Definition 3.12. By the weighted eompaetified amoeba of an algebraic hypersurface H = 
{x GC C** : f{x) ~ 0} we will mean the set pf{H). We denote it by >VC^(/). 


By [9| there exists a deformation p[x) of the Laurent polynomial p{x) such that p{x) — 
p{x) = 8 E M_|_ and moreover the complement of 'WCA{p{x) — e) contains the connected 
component M 3 0 E M** of order 0 G for any e > 0. In other words, there exists 
a deformation of the constant term of p{x) which makes the complement component of 
order zero vanish at the origin. This is illustrated by Fig. ^(d) and 

It is therefore sufficient to check that the origin cannot belong to WCA(p(x)). The 
preimage of the origin under the map pp is given by 


n 


x 


gC 


dp{Q 


dQ 


Cl = \xi\,---,Cn = \Xn\ 


= 0 for any j 
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Since the only point in where the gradient of p{x) vanishes is p, it follows, that 

n = {,..., ) : Xj G [0, 2 %]} . 

To prove the existence of the complement component of order zero we need to show that 
the Laurent polynomial p{x) cannot vanish on 11. 

Consider the set C of lines in such that the intersection of any element in C with S 
contains three points. If jC is empty then the polynomial p{x) is optimal by [7j. From 
now on we assume that C is nonempty. Let L ^ C and consider the restricted polynomial 
p\l{x) := ^ p{s)x^. Since the coefficients of p{x) are defined as the restriction of a 

logarithmically concave function to the support of p(x), the same is true for the coefficients 
of the restricted polynomial p\l{x). Making, if necessary, a monomial change of variables, 
we conclude that p|l(^) = aLm{xY + bLm{x) + c^, where aL^bi^CL G Z and m{x) is 
a monomial function in the coordinates of x. Moreover, this second-order polynomial 
has a logarithmically concave hypergeometric coefficient and vanishes on a certain circle 
centered at the origin. 

The logarithmic concavity of the coefficient of p\l{x) is equivalent to 6^ > XolCl- Since 
the coefficients of p\l{x) are positive integers, its roots are on the same circle centered 
at the origin if and only if p\l{x) = ai^m^xY — 2aLm{p) cosu)m{x) + aLm{pY, where 
coscu < Such a polynomial can only be a restriction of a multivariate hypergeometric 
polynomial in the case when uj = n and its roots are the same. If this holds for any L £ C 
then the initial polynomial p{x) cannot be strongly irreducible. 

This means that there is a connected component of order a in the complement of the 
amoeba Ap(x)- Since a E S C Z” was arbitrary, the amoeba is optimal. □ 

We stress once again that by Definition |3.7| the optimality of a polynomial is the prop¬ 
erty of its zero locus and not the polynomial itself. For instance, the bivariate hyper¬ 
geometric polynomial (1 — x — is optimal since its zeros form the optimal algebraic 
hypersurface {x -\- y = 1}. 

The class of dense optimal multivariate polynomials with Z”-convex supports is of 
course much wider than the class of hypergeometric polynomials. By HU the coefficient 
of a polynomial only has to be ” logarithmically concave enough” for the polynomial itself 
to be optimal. 

Example (a) shows that the condition of (strong) irreducibility is in general not 
necessary for a hypergeometric polynomial to be optimal. On the other hand, the poly¬ 
nomial in Example |3.9| (c) fails the strong irreducibility condition and is not optimal due 
to the fact that it factors into the product of two different Puiseux polynomials with the 
same amoebas. Throughout the rest of the paper we will without loss of generality only 
consider strongly irreducible hypergeometric polynomials. 

Example 3.13. The bivariate Ore-Sato coefficient 


(p(s, t) = (r(t -|- l)r(l + 6 s — 3t)r(31 — 6 s — 2t)) ^ 
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defines a confluent holonomic hypergeometric system with the polynomial solution 

Po{x, y) = l + 593775T + 86493225^^ + 86493225x3 + 593775x4 + x^ + 39331656000x1/+ 
34936343442000x2// + 55898149507200x3// + 216324108000x4// + 54513675216000x//2+ 
2112950051372160000 x 2//2 + 6867087666959520000x3//2 + 10357598291040000x4//2+ 
15382276373989324800000 x 2//3 + 169205040113882572800000x3//3+ 
33807200821954560000 x 4//3 + 3045690722049886310400000x2//4+ 
639595051630476125184000000 x3//4 + 184203374869577124052992000000x3//3+ 
368406749739154248105984000000x3//6. 

The support of this hypergeometric polynomial is Z^-convex and has a triangular convex 
hull. The amoeba of po{x, y) is optimal, see Fig. 




Figure 5. The Newton polygon and the amoeba of po{x,y) 

Remark 3.14. Recall that the Bergman kernel of a complex ellipsoidal domain is given 
by a rational hypergeometric function [13]. One can check that the numerators of such 
rational functions are not necessarily optimal polynomials. The amoebas of the singular 
divisors of the GKZ-hypergeometric functions | 2 ] are known to be solid [T3|. Thus the 
optimal property of the divisors of hypergeometric polynomials cannot be extended to 
the classes of rational or algebraic hypergeometric functions. 

Recall that the Hadamard power of order r G M of a polynomial /(x) = E QjgX 

ses 

is defined to be /^’’^(x) E OgX^. We observe that the set-theoretical limit S{f) := 

ses 

lim C Aff is an amoeba-like simplicial complex. This simplicial complex 

r—)-oo 

for the Hirzebruch polynomial is depicted in Fig. ^ (c) inside the Newton polygon of 
that polynomial. An approximation of the simplicial complex S{pz{x^y)) is depicted in 
Fig. 1^ (c). The geometry of S{f) is related to the amoeba of / while the combinatorics 
of S{f) reflects intrinsic algebraic properties of this polynomial. 

4. Classical bivariate hypergeometric polynomials 

Despite varying terminology, the classical hypergeometric series Fi,..., F 4 , Ci,..., G 3 , 
Hi,... ,H 7 as well as other entries of the Horn list [ 6 ] are universally considered to be 
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intrinsically hypergeometric. For resonant parameters [18], many of these series terminate 
and tnrn ont to be bivariate hypergeometric polynomials. 

Appell’s Fi is one of the most important classical hypergeometric series since by the 
results of | 6 | any bivariate hypergeometric system of second-order equations and holonomic 
rank 3 can be transformed into the system for Fi or a particular limiting case of this 
system. The following statement follows from Theorem |3.10| 

Corollary 4.1. The polynomial instances of the Appell Fi{a, 6 i, 62 , c; x, y) hypergeometric 
function are optimal for a, 61 , 62 , —c < 0 and a > bi + 62 - 

Proof. The imposed conditions on the parameters of Fi(a, 61 , 62 , c; x, y) yield a one-to- 
one correspondence between the F-factors in the coefficient of the power series expansion 
of Fi and the sides of the Newton polygon of its polynomial instance in question. This 
polynomial is therefore hypergeometric in the sense of Definition |3.7| □ 

In Fig.|^we depict the amoeba of the optimal hypergeometric polynomial Fi(—5, —4, —4, 
3;x,?/). 




Figure 6. The Newton polygon and the amoeba of the Appell polynomial 
Fi(-5,-4,-4,3;x,?/) 

Observe however that not every polynomial instance of Fi(a, 61 , 62 , c; x, y) is optimal. 
The optimal property is in general not possessed by the Fi polynomials whose Newton 
polytopes do not have sides that are orthogonal to the gradients of the linear forms in 
the defining Ore-Sato coefficient. For instance, Fi(—4, 5, — 7 , 9 ; x, r/) is not an optimal 
polynomial, its Newton polygon being just a triangle. 

We further remark that the zero locus of a rational instance of a classical hypergeometric 
function need not be an optimal hypersurface. For example, the numerator of the rational 
function 7 ^ 2 ( 5 ; 3/2,1; —1/2, 2; x,y) is not an optimal polynomial. 

5. EXAMPLES 

In this section we collect examples of multivariate hypergeometric polynomials together 
with their Newton polytopes and amoebas. 

Example 5.1. The hypergeometric Horn system defined by the Ore-Sato coefficient 

(p{s, t) = r(s — 6)r(s + t — io)r(t — 6)r(—s -i- 1 — 4)r(— s)r(—s —1 + 2)r(—t)r(s — t — a) 
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admits the following polynomial solution: 

Pi {x, y) = 21x^ + 64a:^ + 21a:^ + 12&xy + 2016a:^y + 4704x^y + 2016 x^ 1 / + 126 x^ 1 / + 21y‘^ + 
2016x1/2 + 22050 x 2|/2 + 47040x3|/2 + 22050x^1/2 + 2016x^?/2 ^ 21x6?/2 + 64i/3 + 4704x?/3 + 
47040x2//^ + 98000x^1/^ + 47040x^//^ + 4704x^i/^ + 64x®i/^ + 21//^ + 2016x//^ + 22050x2//^ + 
47040 x 3//4 + 22050x^1/^ + 2016x^i/^ + 21x®//^ + 126xi/^ + 2016x2//5 + 4704x^1/^ + 2016x^1/5 + 

126x®//^ + 21x2//® _|_ g43;3^6 _|_ 21x^y®. 

(The system itself is too cumbersome to display and we omit it.) The Newton polygon 
and the amoeba of pi{x,y) are shown in Fig. This amoeba turns out to be optimal. 




Figure 7. The Newton polygon and the amoeba of pi{x,y) 


Example 5.2. The next example shows that the number of F-factors in the Ore-Sato 
coefficient of an optimal hypergeometric polynomial can be strictly smaller than the num¬ 
ber of faces of its Newton polytope. The hypergeometric system defined by the Ore-Sato 
coefficient 

(5.1) (p(s, t) = F (s -|- 2t — 5) F (—2s — t — 4) F (—s — 5t -|- 1) 

has the following polynomial solution: 

P 2 {x, y) = 2421619200x®+172972800x®+2882880x^+14560x®+20x9+174356582400x2//+ 
48432384000x3// + 2421619200x^// + 34594560x®// + 160160x®// + 208x^// + 2421619200x//2 + 
691891200x2//2+21621600 x 3//2 + 160160x4//2+286x®//2 + 524160x//3 + 14560x2//3+ hQx^y^ + 
32//^ -|- x//^. 

The support of P 2 {x, y) (bounded by the singular divisors of the corresponding Ore-Sato 
coefficient) and its amoeba are depicted in Fig. 

Example 5.3. The bivariate hypergeometric polynomial supported in the quadrilateral 
with the vertices (2,0), (3, 2), (2,3) and (0,1) is given by pz{x, y) = 240x2 + 3y + 240x// -|- 
1080x2// -|- 30x//2 -|- 180x2//^ -|- 36x^y^ + 2x‘^y^. Fig. (a-c) show the affine amoeba Ap ^, the 
compactified amoeba of P3(x, //), the weighted compactified amoeba of the 6th Hadamard 
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Figure 8. The support of p 2 {x,y) (bounded by the singular divisors 
of (5.1)) and its amoeba 


power of p^{x,y). Fig. shows the vanishing connected component with the order 

(2, 2) in the complement of the weighted compactified amoeba of a deformed version of 
Pz{x,y). We remark that the small component vanishes exactly at the point with the 
coordinates (2, 2), that is, at the order of this component. 



Figure 9. The amoeba, compactified amoeba, the weighted compactified 
amoeba of the 6th Hadamard power and a vanishing connected component 
of the complement to the WCA{p 3 {x,y)) 


Example 5.4. The first maximal minor of the Toeplitz matrix 

/a:i/10000\ 

1 X y 1 0 0 0 

0 1 X 1/ 1 0 0 

0 0 1 X 1/ 1 0 

0 0 0 1 X y 1 

VoOOOlxy/ 

is the degree 6 bivariate Chebyshev polynomial of the second kind [1]. It is optimal in 
the coordinates ^ = xy, p = y^/x which make it dense. 
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